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ABSTRACT 

Estimating  the  Solar  Irradiance  of  an 
intermountain  Region  Using  GOES  Satellite  Data 

by 

Mark  Steven  Walters 
Utah  State  University,  1987 

Maj>  .  ?rofessor:  Dr.  Gail  Bingham 
Department:  Soil  Science  and  Biometeorology 

The  performance  of  two  statistical  models  that  use 
satellite  data  to  calculate  the  global  solar  radiation 
incident  upon  the  earth's  surface  are  assessed.  The 
estimates  are  determined  for  a  mid£latitude  ten  station 
network  and  represent  a  variety  of  sky  cover  conditions. 

Evaluations  of  the  models  for  different  sky  conditions 
reveal  the  need  for  revised  regression  coefficients  for  the 
Hay  and  Hanson  (1978)  model  and  the  Tarpley  (1979)  model. 

The  Hay  and  Hanson  (1978)  model  was  shown  to  perform  better 
for  partly  cloudy  and  overcast  sky  conditions  while  the 
Tarpley  (1979)  model  performed  better  under  clear  skies. 

On  an  hourly  and  daily  time  scale,  the  Hay  and  Hanson  (1978) 
model,  proved  to  be  the  better  performer.  v 


(152  pages) 


CHAPTER  I 


INTRODUCTION 

Surface  solar  radiation  is  of  considerable  significance 
in  fields  as  diverse  as  meteorology,  forestry,  agriculture, 
and  glaciology.  In  addition  to  providing  for  the  basic 
heating  and  cooling  that  generates  the  circulation  of  the 
earth's  atmosphere  and  oceans,  incoming  solar  radiation 
(insolation)  is  responsible  for  the  production  of  green 
plant  foods,  for  providing  an  alternative  energy  source, 
for  activating  the  earth-atmosphere  hydrological  cycle,  and 
for  providing  a  general  environment  suitable  for  human 
habitation.  The  amount  of  insolation  can  be  considered  a 
fundamental  measure  of  the  free  energy  available  at  the 
earth-atmosphere  boundary. 

Current  ground-based  pyranometer  measurements  of 
irradiance  are  limited  to  a  few  weather  stations,  widely 
scattered  universities  and  agricultural  experiment  stations. 
The  quality  of  data  produced  from  this  coarse  network  relies 
on  good  calibration,  regular  maintenance,  and  the  continuous 
functioning  of  all  the  instruments.  These  requirements  are 
not  always  met.  In  addition,  the  data  from  this  network 
can  only  provide  information  about  large-scale  (several 
hundred  kilometers)  variability  and  is  of  little  use  for 
monitoring  smaller  scale  variabilities  or  remote  locations. 


Therefore,  the  current  network  is  generally  insufficient 
to  produce  accurate  insolation  estimates  for  large  areas 
or  remote  locations  of  interest. 
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The  only  practical  sources  of  data  with  the  required 
resolution  and  coverage  are  meteorological  satellites 
such  as  the  Geostationary  Operational  Environmental 
Satellite  (GOES) .  A  GOES  satellite  can  cover  large  areas 
with  adequate  ground  resolution  (1.3  km  at  4 ON)  and  fre¬ 
quency  of  observation  (11-13  Visual  images  per  day)  to  be 
used  for  the  determination  of  solar  irradiance  at  the 
earth's  surface. 

The  possibility  of  using  satellite  data  to  estimate 
solar  irradiance  at  the  surface  has  been  demonstrated  by 
Hay  and  Hanson  (1978),  Tarpley  (1979),  Gautier  et  al . 
(1980),  Brakke  and  Kanemasu  (1981),  Gautier  (1982,  1983), 
Gautier  and  Katsaros  (1984),  Halpern  (1984),  Moser  and 
Raschke  (1984),  Cano  et  al.  (1986),  and  Justus  et  al. 

(1986) . 

The  computer  models  necessary  for  estimating  incoming 
solar  radiation  from  satellite  data  must  consider  the 
interaction  of  radiation  with  the  atmosphere  and  the 
underlying  surface.  This  radiative  transfer  problem  in¬ 
volves  the  absorption,  scattering  and  reflection  of  radia¬ 
tion  energy.  These  processes  determine  the  transfer  of 
radiation  in  the  earth-atmosphere  system  and  are  influenced 
by  the  abundance  of  atmospheric  gases  and  the  concentration 
of  aerosols  and  clouds  as  a  function  of  height. 
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1 .  Objectives 

Currently,  the  irradiance  measurements  for  Utah  are 
obtained  at  26  ground  stations  comprising  the  Utah  Weather 
Network.  The  density  of  surface  measurements  collected 
by  this  network,  while  not  great,  is  still  quite  unusual. 
However,  the  State  Climatologist  for  Utah  has  stated  the 
current  network  fails  to  accurately  quantify  the  solar 
resource  for  over  50%  of  the  state  of  Utah.  The  purpose 
of  this  study  is  to  determine  if  a  reliable  method  of 
estimating  solar  irradiance  using  satellite  data  can  be 
demonstrated  for  this  region. 

The  specific  objectives  of  this  study  are:  1)  To 
assess  the  performance  of  two  published  statistical  models 
that  use  satellite  data  to  estimate  the  global  solar 
irradiance  incident  upon  the  earth's  surface?  and  2)  To 
consider  factors  that  can  be  included  in  these  models  to 
improve  their  use  in  a  mountainous  region  such  as  Utah. 
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CHAPTER  II 

LITERATURE  REVIEW 

The  first  attempts  to  estimate  surface  insolation  that 
involved  the  use  of  satellite  imagery,  utilized  a  variety 
of  satellite  systems.  The  first  satellites,  used  specifi¬ 
cally  in  a  study  of  the  earth's  radiation  budget  and  for 
insolation  estimates,  were  a  series  of  polar  orbiting 
meteorological  satellites. 

The  use  of  polar  orbiting  satellites,  however,  pre¬ 
sented  several  limitations  to  these  early  insolation 
studies.  This  type  of  satellite  can  provide  only  one 
visible  image  per  day  of  any  given  site.  The  hourly  varia¬ 
tion  in  cloudiness,  which  is  the  main  contributer  to  the 
absorption  and  reflectance  of  solar  energy,  could  not  be 
addressed.  Instead,  the  early  work  focused  on  annual, 
seasonal  or  monthly  variations  over  latitudinal  or  zonal 
regions . 

The  question  of  cloud  cover  variability  from  hour  to 
hour  can  be  attempted  with  geostationary  satellite  data. 
This  type  of  data  can  provide  relatively  high  temporal  and 
spatial  resolutions.  This  makes  it  ideal  for  estimating 
incident  solar  radiation  at  the  surface. 

The  first  major  efforts  to  estimate  incident  solar 
radiation  using  geostationary  satellite  data  is  represented 
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by  the  work  of  Hay  and  Hanson  (1978) ,  Tarpley  (1979) , 

Brakke  and  Kanemasu  (1981),  Hiser  and  Senn  (1980),  and 
Gautier  et  al.  (1980). 

The  computer  models  necessary  for  estimating  incident 
solar  radiation  at  the  surface  using  satellite  data  can 
be  divided  into  two  categories:  a  statistical  approach 
based  on  fitted  relationships  between  satellite  and  ground 
measurements,  and  a  physical  approach  using  radiative 
transfer  models  to  formulate  the  relationship  between  the 
satellite  and  ground  measurements. 

Statistical  models  are  more  likely  to  be  precise  when 
applied  to  small  areas,  are  usually  not  as  general,  and 
require  comparison  with  ground  data.  The  work  of  Tarpley 
(1979)  is  perhaps  the  best  known  using  this  method.  The 
physical  approach,  on  the  other  hand,  requires  a  simplified 
model  because  satellites  can  only  measure  a  few  parameters 
among  the  many  that  affect  solar  irradiance.  The  investi¬ 
gations  done  by  Gautier  et  al.  (1980)  are  probably  the 
best  known  of  the  physical  methods. 

On  the  few  models  in  use  today,  most  have  adopted  the 
statistical  approach.  For  example,  Tarpley  (1979),  uti¬ 
lized  a  regression  technique  that  related  satellite  pixel 
brightness  to  insolation.  This  model  estimated  hourly  and 
daily  summer  insolation  values  over  the  Great  Plains  of 
the  United  States.  The  standard  error  of  the  satellite- 
derived  daily  insolation  values  when  compared  against 
pyranometer  values  was  10%  of  the  mean. 


Brakke  and  Kanemasu  (1981)  followed  a  similar  approach 
to  predict  insolation  over  Texas  (Winter  1976)  and  the 
Great  Plains  (Summer  1977) .  Their  technique  produced 
differing  results  for  the  two  seasons.  The  results  for 
the  Winter  1976  data  set  were  within  36%  of  the  observed 
mean,  while  Summer  predictions  were  within  11%  of  the 
Summer  observed  mean. 

Hay  and  Hanson  (1978)  used  hourly  solar  irradiance 
data  from  four  ships,  located  in  the  Atlantic  Ocean,  to 
develop  a  statistical  relationship  between  the  visible 
radiance  data  from  the  GOES  platform  and  the  transmissivity 
of  the  atmosphere  to  solar  radiation.  Consistent  results 
were  found  and  independent  verification  showed  that  six 
hour  irradiance  values  could  be  calculated  to  within  15% 
of  measured  data,  decreasing  to  8%  for  daily  estimates. 

A  revised  version  of  this  model  was  later  used  to  determine 
hourly  and  daily  solar  irradiance  values  over  southwest 
Canada  (Raphael  and  Hay,  1984)  with  similar  results. 

Gautier  et  al .  (1980)  developed  a  simple  physical 
model  taking  into  account  the  effects  of  Rayleigh  scatter¬ 
ing  and  water  vapor  absorption,  but  the  main  emphasis  was 
on  cloud  effects.  Cloud  albedo  and  absorption  were  derived 
from  pixel  data  on  the  assumption  that  they  are  linearly 
related  to  the  cloud  brightness.  Comparisons  with  daily 
insolation  measurements  from  three  pyranometers  located 
in  southern  Canada  showed  that  the  satellite-derived  esti¬ 
mates  were,  on  the  average,  within  9%  of  the  spring  and 
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summer  ground  measurements  for  a  large  variety  of  cloud 
conditions.  The  hourly  variations  monitored  by  the  satel¬ 
lite  also  followed  very  closely  the  variations  measured  on 
the  ground. 

The  above  models  were  developed  and  tested  under  a 
variety  of  radiation  conditions  and  in  differing  environ¬ 
ments:  Tarpley  (1979)  and  Brakke  and  Kanemasu  (1981),  the 

U.S.  Great  Plains;  Hiser  and  Senn  (1980)  utilized  environ¬ 
mentally  different  ground  data  stations  over  the  United 
States;  Hay  and  Hanson  (1978),  the  Tropical  Atlantic;  and 
Gautier  et  al .  (1980),  South  Central  Canada.  Few,  if  any, 
have  attempted  to  apply  their  techniques  to  mountainous 
terrain  where  the  spatial  variation  in  surface  insolation 
can  be  extremely  large. 

The  work  of  Justus  et  al.  (1986),  Klink  and  Dellhopf 
(1986),  Pinker  and  Cario  (1984),  Halpern  (1984),  Cano  et  al. 
(1986),  Moser  and  Raschke  (1984),  and  Powell  et  al .  (1984), 
represent  the  recent  efforts  in  solar  radiation  studies,  all 
using  geostationary  satellite  data. 

These  most  recent  efforts  have  progressed  using, 
generally,  the  same  techniques  as  used  in  the  early  studies. 

A  statistical  method  was  used  in  the  majority  of  the 
studies.  However,  the  study  area  and  the  amount  of  area 
covered  in  each  study  seems  to  vary  widely.  Halpern  (1984) 
developed  a  physical  model  whose  estimates  were  compared  to 
surface  observed  data  from  Northern  California.  Justus 
et  al.  (1986)  used  a  statistical  approach  to  estimate  daily 


total  insolation  on  a  horizontal  surface  at  1“  spacings  in 
latitude  and  longitude  for  the  continental  United  States, 
Mexico  and  parts  of  South  America.  Cano  et  al.  (1986) 
presented  a  statistical  method  for  the  determination  of 
the  global  solar  radiation  of  Europe  using  geostationary 
METEOSAT  images.  Moser  and  Raschke  (1984)  computed  daily 
sums  of  the  downward  solar  radiation  from  METEOSAT  I  and 
II  imagery  over  Europe  and  North  Africa. 


CHAPTER  III 


DATA  ANALYSIS 


1 .  Stratec 


The  strategy  developed  for  this  study  involved  collect 
ing  a  surface  based  ground  truth  data  set  and  a  GOES  satel 
lite  data  set,  for  the  months  of  June  and  August,  1986. 

The  surface  data  set  was  obtained  at  Utah  State  University 
while  the  satellite  data  set  was  obtained  from  the  Atmo¬ 
spheric  Science  Center  at  Colorado  State  University.  Once 
the  data  sets  were  obtained,  two  published  statistical 
models  were  tested  by  comparing  the  calculated  irradiance 
values  to  the  ground  truth  surface  data.  As  a  result  of 
this  comparison,  revised  regression  coefficients  were 
developed  using  three  days  of  data  from  four  of  the  study 
sites.  The  days  used  to  develop  the  regression  coef¬ 
ficients  were  selected  to  represent  clear  and  partly 
cloudy  sky  conditions.  The  revised  models  were  again 
tested  over  the  network  of  sites. 


2.  Data  sets  and  data  processinc 


The  summer  time  period  selected  for  this  study  was 
based  mainly  on  the  availability  of  the  meteorological 
and  satellite  data  sets.  The  meteorological  data  set 


presented  the  least  amount  of  concern  since  reliable  data 


has  been  collected  and  archived  here  at  Utah  State  Uni¬ 
versity  since  1981.  By  far,  the  most  limiting  factor 
presented  during  the  data  collection  phase  of  this  study 
was  the  availability  and  cost  of  the  satellite  data. 

Satellite  data  are  archived  on  a  limited  basis  at  the 
Atmospheric  Science  Center  at  Colorado  State  University. 

A  study  called  Project  FIRE  is  currently  under  way  and  is 
designed  to  better  understand  the  characteristics  of  cirrus 
and  .  tratocumulus  cloud  systems.  Project  FIRE  involves 
collecting  GOES  satellite  data  on  a  6  day  on  -  9  day  off 
schedule  for  a  four-year  period  that  began  in  January, 

1986.  Data  collection  was  intensified  during  the  months  of 
June  through  August,  1986  such  that  hourly  images  ranging 
from  sunrise  to  sunset  were  available.  This  study  utilizes 
the  Tune  and  August  1986,  data  set. 

The  data  set  met  the  basic  conditions  necessary  to 
complete  this  study.  The  primary  requirement  called  for  at 
lea:'*-  five  usable  satellite  images  per  day  so  that  an  accu¬ 
rate  estimate  of  daily  solar  irradiance  could  be 
accomplished  (Justus  et  al.,  1986).  The  second  condition 
called  for  a  spatial  resolution  better  than  8  km.  Addition 
all,  1  he  data  set  offered  at  least  18  days  of  hourly  satel¬ 
lite  imagery,  under  differing  sky  conditions,  for  analyses. 
Lastly,  this  data  set  provided  one  very  important  advantage 
it  v.  -  available  on  a  data  exchange  basis. 


3.  Surface  irradiance  data  set 


The  radiation  data  set  for  this  study  were  collected 
from  10  stations  of  a  26  station  network  comprising  the 
Utah  Weather  Network.  The  area  spans  a  550  km  by  160  km 
{9,000  sq.  km)  section  of  Utah  paralleling  the  Wasatch 
and  Uinta  mountain  ranges.  The  locations  of  the  ten 
stations  of  interest  are  shown  in  Figure  1  and  listed  on 
Table  1.  The  network  was  established  in  1981  by  Dr.  Inge 
Dirmhirn  and  Dr.  Leanard  Hall,  and  is  designed  to  provide 
quality  solar  data  (Bingham  et  al.,  1984).  Among  the 
quantities  measured  at  each  site  are  global  solar  irradi¬ 
ance,  air  temperature,  relative  humidity,  wind  speed,  wind 
direction,  maximum  wind  gust,  precipitation,  rain  gauge 
temperature,  and  soil  temperature.  For  this  study  only 
the  temperature,  relative  humidity,  and  solar  radiation 
.lata  are  of  interest. 

At  each  measurement  site,  the  solar  irradiance  is 
measured  by  a  Licor  LI-200SCZ  silicon  pyranometer.  The 
spectral  response  of  the  silicon  photodiode  extends  over 
a  spectral  region  of  0.4  -  1.1  urn.  Temperature  values  are 
measured  with  a  Campbell  Scientific  101  temperature  probe 
and  relative  humidity  is  measured  with  a  Campbell  Scien¬ 
tific  201  relative  humidity  sensor.  The  instruments  are 
sampled  every  minute  and  recorded  on  a  Campbell  Scientific 
CR2 1  data  logger  and  are  summed  to  provide  hourly  totals. 
The  general  accuracy  of  the  instrumentation  is  well  within 
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Table  1.  Location  of  Utah  Weather  Network  Stations 
used  in  the  study. 


Location 

Latitude 

(N) 

Longitude 

(W) 

Elevation 

(ft) 

Cedar  City* 

37 

45 ' .00" 

113 

01' .50" 

5807 

Delta 

39 

25  '  .45" 

112 

37'  .30" 

4659 

Garland* 

41 

44  '  .00" 

112 

10' .30" 

4400 

Milford* 

38 

22  '  .42" 

112 

59' .00" 

5000 

Springville* 

40 

07' .40” 

111 

42' .30" 

4800 

St.  George 

37 

04  '  .00" 

113 

30' .55" 

2887 

Santaquin 

39 

59’ .55" 

111 

46' .25" 

4800 

Willard 

41 

26'  .44" 

112 

02' .00" 

4350 

Park  City 

40 

39' .00" 

111 

30* .00" 

6800 

Logan 

41 

45' .20" 

111 

47' .30" 

4888 

*  Used  to  develop  revised  regression  coefficients. 


their  published  specifications.  The  data  for  the  10 
stations  of  this  study  are  collected  nightly,  via  tele¬ 
phone  lines,  quality  checked  and  archived  on  a  floppy  disk. 

Factors  that  contribute  to  the  amount  of  insolation 
received  at  the  surface  are  season  of  the  year,  latitude, 
elevation,  time  of  day,  air  quality  (turbidity,  etc.),  and 
cloudiness.  Since  cloud  cover  controls,  to  a  large  extent, 
the  amount  of  solar  radiation  reaching  the  ground,  days 
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from  each  month  were  selected  representing  clear,  partly 
cloudy  and  overcast  cloud  cover  conditions.  Initially,  a 
data  set  of  27  summer  days  was  selected  for  analysis. 

This,  later,  had  to  be  cut  to  18  days  (see  Table  2)  because 
of  navigational  problems  due  to  the  movement  of  the 
satellite . 

4 .  Satellite  data  set 

The  satellite  data  set  used  in  this  study  were  from  a 
Geostationary  Operational  Environmental  Satellite  (GOES) 
operated  by  the  National  Oceanic  and  Atmospheric  Adminis¬ 
tration  (NOAA) .  The  GOES  satellite  was  maintained  at  Earth 
synchronous  altitude,  35800  km  (21480  miles)  above  the 
Earth's  equatorial  plane.  At  this  altitude  its  west  to 
east  motion  equals  that  of  the  Earth  beneath,  ideally 
remaining  stationary  at  a  desired  longitude  (Clark, 

1983)  . 

The  environmental  sensor  onboard  the  satellite  is  a 
Visible  and  Infrared  Spin  Scan  Radiometer  (VISSR) .  The 
sensor  measures  radiance  reflected  from  the  Earth  in  the 
visible  (0.55  -  0.75  um)  and  infrared  (10.5  -  12.6  um) 
regions  of  the  electromagnetic  spectrum  at  a  ground 
resolution  of  1  km  (visible  data)  at  the  equator.  For  this 
study  only  the  visible  data  were  used. 

The  visible  data  are  in  the  form  of  8-bit  count  values 
ranging  from  0-255  counts.  The  minimum  value  is  the  signal 


Table  2.  The  days  and  sky  conditions  used  in  this  study. 


Day/Month 

Julian  Day 

Sky  Condition 

5  June 

156 

partly  cloudy 

6  June 

157 

clear 

7  June 

158 

clear 

8  June 

159 

partly  cloudy 

9  June 

160 

partly  cloudy 

21  June 

172* 

clear 

22  June 

173 

clear 

23  June 

174* 

partly  cloudy 

24  June 

175 

partly  cloudy 

25  June 

176 

partly  cloudy 

5  August 

217 

partly  cloudy 

6  August 

218 

clear 

8  August 

220 

partly  cloudy 

9  August 

221 

clear 

10  August 

222 

clear 

20  August 

232 

partly  cloudy 

21  August 

233* 

partly  cloudy 

24  August 

236 

partly  cloudy 

*  Used  to  develop 

revised  regression 

coefficients . 

received  from  a  black  surface,  and  the  maximum  value  is  the 
signal  received  from  a  surface  of  100%  reflectance.  The 


surface  is  also  assumed  to  reflect  radiation  incident  on 
it  equally  in  all  directions.  This  type  of  surface  is 
called  a  Lambertian  surface. 

The  satellite  data  used  in  this  study  are  archived 
at  the  Atmospheric  Science  Center  of  the  University  of 
Colorado.  The  data  are  stored  in  digital  form  on  magnetic 
tapes  and  can  be  used  both  as  image  and  digital  data. 
Filtered  satellite  data  corresponding  to  the  days  listed 
in  Table  2  were  obtained  by  first  specifying  the  number 
of  pixel  lines  and  elements  in  the  image  (in  this  case 
512  by  512),  and  by  specifying  the  latitude  and  longitude 
coordinates  of  the  image  center.  All  of  the  images  were 
centered  on  40N,  HOW  and  contain  all  the  stations  in  the 
Utah  Weather  Network  plus  areas  of  Nevada,  Colorado, 
Arizona,  and  Wyoming. 

5 .  Problem  areas 

To  properly  test  the  models  used  in  this  study  the 
meteorological  and  satellite  data  sets  described  in  previous 
sections  must  be  combined.  These  data  sets  are  vastly 
different  leading  to  a  number  of  problems. 

a.  Completeness  of  data  sets 

As  mentioned  in  previous  sections,  the  Project  FIRE 
satellite  data  were  collected  on  a  6  day  on  -  9  day  off 
schedule.  This  schedule  allows  the  collection  of  only 
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twelve  days  of  imagery,  from  each  month,  for  analyses. 
Additionally,  two  images  for  each  day  were  not  collected 
during  Project  FIRE.  They  are  for  the  hours  of  7  AM  and 
6  PM  local  time  (DST) .  For  a  data  perfect  condition, 
twelve  satellite  images  were  available  from  6  AM  to  9  PM 
local  time  (DST).  For  the  month  of  August,  only  ten 
visible  images  were  available  due  to  the  later  sunrise  and 
an  earlier  sunset. 

The  surface  data  for  the  ten  sites  in  this  study  are 
automatically  collected  via  telephone  lines.  This  col¬ 
lection  system  was  installed  during  the  summer  of  1986. 
Seven  of  the  sites  were  on  line  by  June  and  three  came  on 
line  during  July  (Logan,  Park  City,  and  Santaquin) .  June 
data  for  these  three  sites  were  not  available.  Addition¬ 
ally,  equipment  failures  caused  missing  data  at  scattered 
sites  and  times. 

As  a  result  of  missing  data  from  the  satellite  and  the 
surface  irradiance  data  sets,  most  of  the  18  days  evaluated 
in  this  study  have  some  missing  data. 

b.  Spectral  ranges  and 

sensor  viewing  angles 

The  satellite  measurements  were  instantaneous  measure¬ 
ments,  made  once  an  hour,  over  a  small  solid  viewing  angle. 
The  pyranometer  measurements  were  collected  from  all  angles 
of  a  hemisphere  and  were  integrated  over  an  hour.  An 
assumption  is  made  that  an  instantaneous  satellite  measure¬ 
ment  taken  on  the  hour  can  be  made  to  approximate  a 


MV 
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pyranometer  measurement  averaged  over  an  hour.  Addition¬ 
ally,  the  pyranometer  measurements  extend  over  a  spectral 
region  of  0.4  -  1.1  um  while  the  satellite  data  extend 
over  the  smaller  visible  band  (0.55  -  0.75  um) .  An  attempt 
is  made  to  partly  account  for  the  angle  and  time  discrepan¬ 
cies  by  spatially  averaging  the  satellite  radiation  esti¬ 
mates  over  a  5  x  5  pixel  box.  Averaging  over  the  5x5 
array  also  attempts  to  account  for  the  variations  in  the 
eight  visible  sensors  on  the  VISSR  (Gautier,  1982).  The 
5x5  pixel  array  size  was  based  on  a  navigational  accuracy 
of  one  to  two  pixels. 

c.  Navigation 

The  accuracy  of  the  results  obtained  in  this  study  are 
directly  related  to  the  ability  to  align  points  on  a  satel¬ 
lite  image  to  the  same  points  on  the  Earth's  surface.  This 
is  accomplished  by  means  of  a  procedure  called  "navigation." 

To  assure  accurate  alignment  of  the  imagery,  navigation 
was  accomplished  in  two  separate  procedures.  Initial 
navigation  of  the  satellite  imagery  was  performed  using  a 
model  described  by  Hambrick  and  Philips  (1980).  The  model 
is  based  on  the  knowledge  of  some  physical  parameters  of 
the  satellite  such  as  orbit,  altitude,  etc.,  and  both 
stellar  and  terrestrial  navigation  points.  The  relative 
accuracy  of  this  procedure  is  claimed  to  be  one  pixel. 
However,  visual  comparison  of  the  imagery  on  the  COMTAL 
image  processor  often  revealed  much  larger  errors. 


19 


To  accurately  navigate  the  satellite  imagery,  a  final 
procedure  called  "roaming"  was  accomplished.  This  pro¬ 
cedure  relies  totally  on  visual  identification  of  land¬ 
marks  on  the  Earth's  surface  and  depends  on  relatively 
clear  skies.  An  accuracy  of  one  or  two  pixels  was  achieved 
by  selecting  a  particularly  clear  and  sharp  image  as  a 
"base  image,"  then  visually  inspecting  each  satellite  image 
on  the  COMTAL  image  processor  and  aligning  the  positions 
of  known  landmarks  to  those  of  the  base  image.  When  se¬ 
quential  images  were  compared,  a  consistency  of  position 
of  a  given  feature  was  achieved.  Examples  of  landmarks 
used  to  navigate  the  images  extend  over  Utah,  Nevada,  and 
Colorado  and  include  the  Great  Salt  Lake,  Lake  Mead,  Lake 
Powell,  Bear  Lake,  the  Great  Salt  Lake  Desert,  and  the 
Wasatch  and  San  Juan  Mountain  Ranges  (see  Figure  1).  In 
clear  sky  conditions  the  primary  navigational  features 
were  the  Great  Salt  Lake,  Lake  Mead,  and  Lake  Powell.  In 
nearly  every  image  at  least  two  of  the  primary  navigational 
features  were  visible.  It  was  necessary  in  only  a  few  of 
the  overcast  situations  to  use  secondary  navigational 
features  such  as  mountain  ranges  and  smaller  lakes. 

An  additional  and  separate  navigational  problem  was 
discovered  during  the  satellite  data  collection  phase. 

The  current  GOES  network  consists  of  only  one  operational 
GOES  satellite,  instead  of  two.  The  GOES-East  (GOES-5) 

VT.SSR  ceased  to  operate  in  the  Spring  of  1985.  As  a  result, 
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the  only  operational  satellite  (GOES-6)  was  moved  during 
the  hurricane  season  of  1986  to  better  cover  the  western 
Atlantic.  It  took  approximately  40  days  to  move  the 
satellite  from  its  primary  location  at  108W  to  its  summer 
location  at  98W.  The  satellite  was  moving  from  18  June 
to  28  July  and  created  several  unforeseen  problems  as 
related  to  this  study. 

The  first  problem  concerns  the  accuracy  of  the  navi¬ 
gation  during  the  time  the  satellite  was  in  motion;  from 
June  18  to  28  July.  To  save  fuel,  the  satellite  is  moved 
slowly  so  the  change  in  the  viewing  angle  of  the  satellite 
was  negligible  for  the  first  few  days.  For  this  reason, 
the  images  through  June  24  proved  to  be  acceptable.  The 
data  for  July,  however,  proved  to  be  unusable  due  to  errors 
in  navigating  the  images  as  the  satellite  was  relocated. 

Once  the  satellite  was  at  its  summer  location  of  98W, 
it  became  necessary  to  roam  the  August  satellite  images 
to  a  base  image  selected  from  the  August  imagery.  Clearly, 
once  the  satellite  had  been  relocated,  the  location  of 
every  navigational  feature  would  change  as  well  as  the 
station  pixel  coordinates  (see  Table  3) . 

It  was  decided  to  treat  the  June  and  August  data  as 
separate  data  sets  and  account  for  the  location  of  the 
satellite  within  the  two  models.  Namely,  the  location 
of  the  satellite  for  the  June  data  would  be  108W  and  the 
location  of  the  satellite  for  the  August  data  would  be 


98W. 


Table  3.  The  latitude,  longitude  and  pixel  coordinates 
for  each  measurement  site  during  June  and  August. 


Latitude 

Longitude 

Pixel  Coordinates 

(N) 

(W) 

June 

August 

X 

Y 

X 

Y 

Cedar  City 

37.75 

113.025 

171 

365 

174 

341 

Delta 

39.42 

112.62 

191 

234 

195 

218 

Garland 

41.73 

112.17 

213 

52 

218 

49 

Mi.  Iford 

38.37 

112.88 

173 

317 

177 

296 

Gpringville 

40.12 

111.705 

236 

179 

241 

167 

St.  George 

37.07 

113.51 

147 

419 

150 

391 

Santaquin 

39.99 

111.77 

233 

189 

237 

177 

Willard 

41.44 

112.03 

220 

75 

224 

70 

Park  City 

40.65 

111.50 

246 

137 

250 

128 

Logan 

41.75 

111.78 

232 

50 

236 

47 

The  relocation  of  the  satellite  resulted  in  having  to 
develop  separate  models  for  both  months,  selecting  a 
separate  base  image  for  both  months,  specifying  the  sepa¬ 
rate  station  pixel  locations  for  each  month  and,  of  course, 
Khe  deletion  of  all  the  July  data. 
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6 .  Precipitable  water 


Moisture  is  one  of  the  meteorological  parameters  that 
affects  the  amount  of  solar  radiation  reaching  the  surface. 
The  Tarpley  (1979)  model  accounts  for  this  variable  by 
including  precipitable  water  within  its  framework.  Tarpley 
used  radiosonde  based  precipitable  water  data  from  the 
National  Meteorological  Center  (NMC)  data  file.  The  files 
were  updated  twice  daily;  the  0000  GMT  file  contained 
information  most  nearly  time-coincident  with  the  satellite 
data.  Additionally,  precipitable  water  values  for  each 
target  were  retrieved  from  the  nearest  NMC  grid  point 
which  could  be  as  far  as  2  degrees  latitude  and  longitude 
from  the  target  (Tarpley,  1979).  Obviously,  this  technique 
could  be  open  to  much  error. 

An  alternative  technique  was  used  in  this  study.  Using 
Smith's  (1966)  empirical  formulation,  given  a  surface 
dewpoint  temperature,  precipitable  water  values  can  be 
calculated  for  each  of  the  measurement  stations,  based  on 
a  model  atmosphere.  The  use  of  Smith's  emperical  technique 
offers  the  ability  to  update  the  precipitable  water  data  on 
an  hourly  basis  for  all  of  the  study  sites.  The  published 
equation  is: 

In  (u)  =  [0.1133  -  In ( X  +  1)1  +  0.0393  Td  (Smith,  1966) 

u  =  precipitable  water  in  cm 

Td  =  surface  dewpoint  temperature  (F) 

X  =  latitude  and  seasonally  adjusted  coefficient 


Atwater  and  Ball  (1976)  computed  solar  radiation  for 
eleven  stations  in  the  U.S.  They  showed  that  the  method 
of  determining  precipitable  water,  from  radiosonde  obser¬ 
vations  or  estimating  it  from  the  surface  dewpoint  temper¬ 
ature,  had  little  effect  on  the  calculated  solar  radiation. 
Raphael  (1983)  tested  this  finding  by  calculating  solar 
radiation  using  Smith's  formulation  and  found  precipitable 
water  values  determined  from  the  surface  dewpoint  temper¬ 
ature  produced  solar  radiation  estimates  closer  to  the 
observed  solar  radiation,  and  that  precipitable  water 
values  derived  from  the  radiosonde  data  were  largely 
inappropriate  due  to  the  errors  introduced. 

7 .  The  models 

The  models  chosen  for  evaluation  are  those  of  Hay  and 
Hanson  (1978)  and  Tarpley  (1979).  The  two  models  were 
originally  developed  for  different  environments  and  tested 
under  a  variety  of  radiation  conditions,  but  neither  has 
been  tested  for  use  over  a  mountainous  region  such  as  Utah. 
For  this  reason,  both  models  were  initially  tested  using 
the  original  model  parameters  and  regression  coefficients 
over  select  days  of  the  study  period. 

The  model  developed  by  Hay  and  Hanson  (1978)  is  the 

simplest  of  the  two  models.  The  basic  equation,  where  the 

irradiance  at  the  surface  I  is  written  as: 

s ' 
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I  =  I  cose(a-bSR)  (Raphael,  1983) 

s  o 

I  =  solar  constant  (1353  W/m) 

SR  =  normalized  satellite  reflectance 
a,b=  regression  coefficients 
0  =  local  solar  zenith  angle 

The  Hay  and  Hanson  (1978)  model  utilizes  a  satellite 
calibration  technique  that  converts  the  satellite  pixel 
brightness  counts  into  a  reflectance  value.  The  calibra¬ 
tion  procedure  used  in  the  model  was  provided  by  E.  Smith, 
Colorado  State  University  (Raphael,  1983).  Using  this 
procedure  a  computer-based  "look-up  table"  is  generated 
relating  the  brightness  counts  to  relevant  normalized 
reflectance  and  irradiance  values.  The  normalized 
reflectance  for  any  pixel  brightness  value  is  obtained 
simply  by  locating  this  brightness  value  in  the  look-up 
table  and  retrieving  the  appropriate  normalized  reflectance 
(Raphael  and  Hay,  1984) . 

The  value  retrieved  from  the  look-up  table  is  used  in 
the  final  form  of  the  irradiance  calculation.  The  final 
form  of  the  irradiance  calculation  is  written  as: 

I  =  a(xext)  +  b(xext)xir 
a,b  =  regression  coefficients 

xext  =  the  extraterrestrial  global  irradiance  that 
is  calculated  from  the  product  of  the  solar 
constant  and  the  cosine  of  the  zenith  angle 
at  the  satellite  image  time. 
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xir  =  the  product  of  the  quantity  retrieved  from 

the  look-up  table,  the  normalized  reflectance 
value,  and  the  inverse  of  the  cosine  of  the 
solar  zenith  angle. 

The  only  supplementary  data,  the  surface  measured 
irradiance  value,  is  used  by  the  model  in  a  statistical 
comparison  to  the  calculated  irradiance  value. 

The  model  presented  by  Tarpley  (1979) ,  was  developed 
and  tested  using  data  from  the  United  States  Great  Plains. 
An  important  feature  of  the  model  is  the  brightness  param¬ 
eterization  given  by: 

2 

B  =  a  +  bcose  +  csin0  cos<j>  +  dsin0  cos  <|>  (Tarpley, 

B  =  predicted  minimum  brightness  1979) 

<J>  =  azimuth  angle  between  sun  and  satellite 
6  =  local  solar  zenith  angle 
a,b,c,d  =  regression  coefficients 
This  equation  accounts  for  the  changing  incident  flux, 
the  changes  in  the  target  brightness  due  to  shadowing  at 
the  surface  and  anistropic  scattering. 

Cloud  coverage  is  determined  from  a  two-threshold 
method  as  presented  by  Shenk  and  Salomonson  (1972)  .  Three 
categories  are  determined  by  this  method.  They  are  clear, 
partly  cloudy  (50%  cloud  cover)  and  cloudy  (100%  cloud 
cover) .  The  cloud  factor  (n)  is  computed  using  the 
equation : 


•5N2  +  N3 


(Tarpley,  1979) 


N2  +  2N3 

Nx  +  N2  +  N3  2N 

n  =  cloud  factor 
N  =  total  number  of  pixels  in  target  area 
1^,  N2,  N3  =  number  of  pixels  in  clear,  partly  cloudy, 
and  cloudy  categories  respectively 
■’he  clear/partly  cloudy  threshold  (Tl)  is  the  predicted 
cl  r  brightness  plus  twelve  counts.  Any  pixel  value 
greater  than  or  equal  to  Tl  is  considered  clear.  The 
paialy  cloudy/cloudy  threshold  (T2)  is  the  predicted  clear 
brightness  plus  twenty  counts.  Any  pixel  value  greater 
th  m  Tl  but  less  than  T2  is  considered  to  be  partly  cloudy. 
Pi  a  :  i  s  greater  than  T2  are  considered  cloudy. 

iiree  regression  equations  were  developed  to  estimate 
tl  i rradiance  at  the  surface  under  clear,  partly  cloudy 
an1  overcast  sky  conditions.  For  clear  sky  conditions,  when 
th.*  loud  factor  (n)  is  less  than  .4,  the  hourly  surface 
iricidiance  ( I s )  is  calculated  using 

2 

tr.  =  ai  +  cose  +  c^  V  +  d^n  +  e^  (Im/B)  (Tarpley, 

1979) 

Wh*  n  a  partly  cloudy  condition  occurs,  the  cloud  factor 
(n)  is  between  .4  and  1.0,  the  hourly  surface  irradiance 
Ic  is  calculated  using 

2 

I  .  =  a„  +  b_  cose +  c0n  (Icld/B  )  (Tarpley,  1979) 

li  z.  z  z  o 

Wh the  cloud  factor  (n)  equals  1.0  an  overcast  sky 
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condition  exists  and  I  is  computed  using  the  equation 

Ig  =  a^  +  b^  cosG+c^  (Icld/BQ)2  (Tarpley,  1979) 

where  I  is  the  hourly  surface  irradiance,  I  is  the  mean 

m 

target  brightness,  B  is  the  predicted  clear  brightness, 
Icld  is  the  mean  cloud  brightness,  Bq  is  the  normalized 
clear  brightness,  the  regression  coefficients  are 
a,b,c,d,e  (the  original  values  are  used  here)  and  ( ^  )  is 
the  atmospheric  transmittance. 

Atmospheric  transmittance  ( ^  )  is  calculated  using  the 
following  equation 

ip  =  ip  •  4>  ‘*P  (Tarpley,  1979) 

ws  wa  r  fit 

=  l-0.00225um  transmission  due  to  water  vapor 
scattering 

r  =  l-0.077(wm)‘  transmission  due  to  water  vapor 
wa 

absorption 

ip  2 

r  =  0 . 972-0 . 0826m  +  0.00933m  transmission  due  to 
Rayleigh  scattering 
u  =  precipitable  water  (cm) 
z  =  station  elevation  (m) 

m  =  optical  air  mass  =  e  / [cosG  +.15/ 

1  25  3 

(93.885-8  ) 


Satellite  azimuth  angles 


£ 


To  determine  the  satellite  azimuth  angles  (a) ,  necessary 
in  both  model  calculations,  several  equations  from  Sellers 
(1965)  are  used.  The  satellite  azimuth  angle  from  south  is 


v.  *y*' 


TTl'V-l' 


calculated  by  first  determining  the  sun's  azimuth  (Z). 
Using  spherical  trigonometry  it  follows: 
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cosZ  =  sin<J>  sin6  +  cost))  cos6  cosh  (Sellers,  1965) 

<}>  =  station  latitude 

6  =  angular  distance  of  the  satellite  north  or  south 

of  the  equator.  For  a  geostationary  satellite  6=  0. 
h  =  hour  angle.  The  longitude  of  the  station  minus  the 
longitude  of  the  satellite. 

The  satellite  azimuth  angle  can  then  be  calculated 
using : 

cos  5  sin  h  , 

sm  a  =  — sTn"Z -  (Sellers,  1965) 

As  mentioned  earlier,  the  movement  of  the  satellite 
from  98W  to  108W  during  the  study  period  required  the 
development  of  two  sets  of  satellite  azimuth  angles. 
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CHAPTER  IV 


RESULTS  AND  DISCUSSION 


1 .  The  Hay  and  Hanson  model 


The  initial  test  of  the  Hay  and  Hanson  (1978)  model  was 
performed  for  all  of  the  study  sites  using  all  the  avail¬ 
able  days  from  the  data  set.  Prior  to  the  test,  Julian 
days  172,  representing  clear  sky  conditions  and  days  174 
and  233  were  selected  for  analysis  representing  partly 
cloudy  sky  cover  conditions. 

The  averaged  hourly  calculations  from  the  initial  test 
indicate  that  the  model  is  unbiased  towards  any  of  the  test 
sites.  The  results  also  indicate  that  the  model  is  gener¬ 
ally  consistent  for  two  sky  cover  situations.  The  model 
underestimates  clear  conditions  that  remain  clear  all  day 
(day  172)  and  overestimates  partly  cloudy  situations  that 
remain  partly  cloudy  all  day  (day  233) .  The  models  per¬ 
formance  seems  to  decrease  for  those  conditions  that  start 
the  day  clear  and  end  the  day  with  cloudy  skies.  The 
values  listed  in  Table  4  and  Figures  2-12  indicate  the 
performance  of  the  model  at  select  locations  of  the  study 
area  (see  Appendix  for  complete  statistics  on  the  individual 
stations) . 


*  v.v,  <Vi\ 


Table  4.  Hourly  statistics,  at  select  sites,  from  the 
initial  test  of  the  Hay  and  Hanson  model  (1978)  using  the 
original  regression  coefficients.  The  three  days  represent 
clear  (day  172),  and  partly  cloudy  (days  174,  233)  sky 
conditions . 


IRRADIANCE 


DAY 

OBS 

CALC 

N 

MBE% 

RMSE% 

LOCATION 

172 

26049.0 

25413.6 

11 

2.5 

8.3 

Cedar  City 

172 

25942.9 

25921.9 

11 

.1 

.3 

St.  George 

172 

26494.1 

26359.5 

11 

.5 

1.7 

Garland 

172 

26305.7 

26283.6 

11 

.1 

.3 

Springville 

174 

20562.5 

22513.0 

11 

-8.7 

28.7 

Cedar  City 

174 

23403.2 

22600.1 

11 

3.6 

11.8 

Milford 

174 

23564.4 

23446.3 

11 

.5 

1.7 

Delta 

233 

12459.4 

14411.3 

9 

-13.5 

40.6 

Milford 

233 

15059.9 

19179.8 

9 

-21.5 

64.4 

Park  City 

233 

14469.1 

17458.3 

9 

-17.1 

51.4 

Logan 

233 

19094.4 

19868.6 

9 

-3.9 

11.7 

Garland 

TOT: 

233404.6 

243456.0 

113 

-4.1 

The  figures 

depicting 

day 

172,  the 

clear 

day,  indicate 

that  the  model  tends  to  underestimate  the  irradiance  during 
mid  day  and  overestimate  during  the  morning  and  evening 
hours.  The  values  in  Table  4  indicate  that  these  high  and 


Calculated  and  observed  irradiance  using  the  original  regres 
coefficients  from  the  Hay  and  Hanson  (1978)  model. 


Fig.  3.  Calculated  and  observed  irradiance  using  the  original  regression 
coefficients  from  the  Hay  and  Hanson  (1978)  model. 


Calculated  and  observed  irradiance  using  the  original  regression 
coefficients  from  the  Hay  and  Hanson  (1978)  model. 


Fig.  5.  Calculated  and  observed  irradiance  using  the  original  regression 
coefficients  from  the  Hay  and  Hanson  (1978)  model. 


Calculated  and  observed  irradiance  using  the  original  regression 
coefficients  from  the  Hay  and  Hanson  (1978)  model. 


Fig.  8.  Calculated  and  observed  irradiance  using  the  original  regression 
coefficients  from  the  Hay  and  Hanson  (1978)  model. 
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Calculated  and  observed  irradiance  using  the  original  regression 
coefficients  from  the  Hay  and  Hanson  (1978)  model. 


Calculated  and  observed  irradiance  using  the  original  regression 
coefficients  from  the  Hay  and  Hanson  (1978)  model. 


Calculated  and  observed  irradiance  using  the  original  regre 
coefficients  from  the  Hay  and  Hanson  (1978)  model. 


Fig.  12.  Calculated  and  observed  irradiance  using  the  original  regression 
coefficients  from  the  Hay  and  Hanson  (1978)  model. 


low  estimates  for  clear  days  don’t  quite  cancel  out, 
leaving  calculated  values,  in  most  cases,  lower  than  the 
observed. 

Figures  9-23,  depicting  partly  cloudy  sky  conditions, 
reveal  that  the  model  generally  mirrors  the  observed  data 
line.  The  figures  also  indicate  that  the  model  tends  to 
overestimate  irradiance  values  during  the  morning  and 
evening  hours  with  underestimation  during  the  mid  day  hours 
producing  a  greater  amount  of  underestimation  during  the 
mid  day  hours  than  with  the  mid  day  under  estimations  for 
clear  days. 

The  figures  depicting  partly  cloudy  days  also  indicate 
that  the  model  matches  the  late  morning,  evening  and  mid 
day  observations  quite  well.  However,  the  model  has 
trouble  with  peaks  and  valleys  in  the  observed  data  often 
missing  major  events  of  cloudiness  or  sunshine.  This  often 
produces  the  most  erroneous  irradiance  calculations  of  the 
three  sky  cover  conditions. 

The  fact  that  the  model  has  a  problem  calculating 
irradiance  values  for  cloudy  days  should  come  as  no  sur¬ 
prise.  It  must  be  remembered  that  the  irradiance  calcula¬ 
tion  is  based  on  one  satellite  image  taken  at  the  beginning 
of  the  hour  while  the  observed  data  is  sampled  every  minute 
and  integrated  for  the  hour.  An  assumption  is  made  that 
the  satellite  measurement  taken  once  an  hour  can  approxi¬ 
mate  a  pyranometer  measurement  averaged  for  the  same  hour. 
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Calculated  and  observed  irradiance  using  the  improved  regression 
coefficients  in  the  Hay  and  Hanson  (1978)  model. 


Calculated  and  observed  irradiance  using  the  improved  regression 
coefficients  in  the  Hay  and  Hanson  (1978)  model. 


Calculated  and  observed  irradiance  using  the  improved  regression 
coefficients  in  the  Hay  and  Hanson  (1978)  model. 


Fig.  16.  Calculated  and  observed  irradiance  using  the  improved  regression 
coefficients  in  the  Hay  and  Hanson  (1978)  model. 


Calculated  and  observed  irradiance  using  the  improved  regression 
coefficients  in  the  Hay  and  Hanson  (1978)  model. 


Fig.  18.  Calculated  and  observed  irradiance  using  the  improved  regression 
coefficients  in  the  Hay  and  Hanson  (1978)  model. 


Fig.  19.  Calculated  and  observed  irradiance  using  the  improved  regression 
coefficients  in  the  Hay  and  Hanson  (1978)  model. 


in  the  Hay  and  Hanson  (1978)  model. 


Fig.  21.  Calculated  and  observed  irradiance  using  the  improved  regression 
coefficients  in  the  Hay  and  Hanson  (1978)  model. 


Hi  (■PUNMIU) 

Gotoh/e-pi  m)  3ohvkjv*wi 

si 

g 


,v  v 


v*. 


■W 


Fig.  23.  Calculated  and  observed  irradiance  using  the  improved  regression 
coefficients  in  the  Hay  and  Hanson  (1978)  model. 
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For  example,  the  instantaneous  measurement  made  by  the 
satellite  may  see  clear  skies  while,  in  fact,  the  pyrano- 
meter  has  been  shadowed  by  drifting  clouds  during  the 
remainder  of  the  hour.  This  would  produce  a  high  calculated 
value  and  a  low  observed  value. 

The  accuracy  of  the  original  model  was  determined  by 
running  a  statistical  analysis  of  data  listed  in  Table  4. 

The  totaled  value  listed  at  the  bottom  of  the  table  indi¬ 
cates  the  model's  performance  under  all  sky  cover  conditions. 
The  original  Hay  and  Hanson  (1978)  model  performed  sur¬ 
prisingly  well  considering  where  the  model  was  developed. 
However,  limitations  in  the  model  become  apparent  when  it 
is  used  to  calculate  cloudy  conditions.  Although  the  per¬ 
formance  of  the  model  is  acceptable  for  clear  days,  an 
attempt  is  made  to  improve  the  models  overall  performance 
by  determining  new  regression  coefficients. 


a.  Rearession  coefficients 


The  development  of  new  regression  coefficients  was 
accomplished  using  data  collected  at  four  of  the  study 
sites.  The  sites  utilized  in  the  analyses  are  Cedar  City, 
Milford,  Garland  and  Springville.  The  sites  were  selected 
to  represent  various  conditions  (elevation,  latitude, 
longitude,  etc.)  over  the  network.  The  analysis  centered 
on  days  172,  174,  and  233  to  represent  various  sky  cover 
conditions.  The  new  coefficients  were  calculated  using  an 


available  software  package  called  Number  Cruncher.  They 
are  determined  to  be :  a  =  0.77  and  b  =  -0.74.  Both 
regression  coefficients  changed  only  slightly  from  the 
values  used  in  the  original  model.  The  old  and  new 
regression  coefficients  are  listed  in  Table  5. 

The  original  Hay  and  Hanson  model  was  developed  for 
use  over  the  tropical  Atlantic  and  predominately  for  use 
during  the  summer  months.  Several  considerations  exist 
when  implementing  the  model  at  mid-latitude  locations  such 
as  Utah.  The  major  considerations  are  latitude,  elevation, 
predominate  cloud  type,  cloud  absorption,  and  lower  amounts 
of  water  vapor.  Of  these  three  considerations,  the  biggest 
contributer  to  variability  in  the  model  could  be  lower 
water  vapor  amounts. 

b.  The  revised  model  results 

The  revised  model  was  also  tested  over  all  the  sites 
using  all  the  available  days  from  the  data  sets.  The  hourly 
average  statistics,  at  select  sites,  for  days  172,  174  and 
233  are  listed  in  Table  6  (see  Appendix  for  complete 
statistics  on  the  individual  stations) . 

The  effect  of  the  revised  regression  coefficients  for 
all  sky  cover  conditions  can  be  seen  when  the  totaled 
irradiance  calculations  from  Tables  4  and  6  are  compared. 

The  mean  bias  error  has  been  considerably  decreased  and 
noticeable  improvement  is  observed  in  most  of  the  model 


calculations . 


Table  5.  The  original  and  revised  regression  coefficients 
as  developed  for  this  study. 


i 

§ 

COEFFICIENT  ORIGINAL  REVISED 


a 

.79 

.77 

ft 

b 

-.71 

-.74 

9 

h  «• 

Table 

6.  The 

hourly  statistics  using 

the  revised  regression 

hit 

coefficients  from  the  Hay 

and 

Hanson 

(1978)  model  for  clear 

£ 

(day 

1 

172),  and 

partly  cloudy 

(days  174,  223)  sky  conditions. 

I? 

IRRADIANCE 

i 

|  DAY 

OBS 

CALC 

N 

MBE% 

RMSE% 

LOCATION 

|  172 

26049.0 

24550.3 

11 

3.8 

13.0 

CEDAR  CITY 

172 

25942.9 

25079.4 

11 

3.4 

11.4 

ST.  GEORGE 

H 

172 

26494.1 

25545.0 

11 

3.7 

12.3 

GARLAND 

ft 

172 

! 

26305.7 

25461.2 

11 

3.3 

11.0 

SPRINGVILLE 

174 

21890.3 

23668.7 

12 

-7.5 

26.0 

CEDAR  CITY 

174 

25847.0 

24559.1 

12 

5.2 

18.2 

MILFORD 

% 

174 

26418.6 

25060.6 

12 

5.4 

18.8 

DELTA 

V 

V 

s. 

233 

12459.4 

13347.9 

9 

-6.7 

20.0 

MILFORD 

ffi 

233 

15059.9 

18346.5 

9 

-17.9 

53.7 

PARK  CITY 

ft 

'  233 

14469.1 

16566.9 

9 

-12.7 

38.0 

LOGAN 

8 

|  233 

19094.4 

19078.8 

9 

.1 

.2 

GARLAND 

TOT:  240030.4  241264.4  116 
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From  the  results  in  Table  6  and  those  listed  in  the 
Appendix,  the  improvement  in  the  MBE%  and  RMSE%  values  for 
partly  cloudy  days  indicate  the  performance  of  the  model 
for  this  type  of  day  has  been  improved.  The  smaller  value 
for  the  coefficients  means  that  the  mid  day  estimations 
would  decrease  with  slightly  higher  estimations  during  the 
morning  and  evening  hours.  This  trend  would  be  most  bene¬ 
ficial  at  sites  that  experience  cloudy  conditions  during 
the  peak  insulation  period. 

Figure  17  shows  such  a  situation.  The  conditions  at 
Cedar  City  on  day  174  indicate  that  the  site  began  the  day 
cloud  free  and  received  the  normal  amount  of  irradiance 
for  such  a  condition.  During  mid  morning,  convective 
clouds  began  to  develop  and  remained  through  the  mid  day 
hours,  not  allowing  the  site  to  receive  the  maximum  amount 
of  irradiance.  The  site  again  became  cloud  free  during 
the  evening  hours  and  then  followed  a  clear  sky  trend.  For 
this  situation,  the  revised  model  produced  estimates  that 
were  slightly  high  during  the  morning  and  evening  hours, 
which  is  similar  to  a  clear  day.  The  mid  day  estimates 
were  similar  to  a  cloudy  day.  The  two  combinations  help 
to  balance  the  estimates  produced  for  this  type  of  day. 

The  days  that  showed  no  improvement  can  be  characterized 
in  two  ways.  First,  the  site  had  widely  divergent  cloud 
conditions  from  hour  to  hour.  In  a  case  such  as  this,  one 
satellite  image  per  hour  would  be  a  poor  indicator  of  the 
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conditions  experienced  at  the  site.  Second,  the  site 
experienced  cloud  free  conditions  during  the  mid  day  hours 
and  was  allowed  to  receive  the  same,  or  nearly  the  same, 
irradiance  as  a  clear  site.  Most  of  the  sites  that  show 
no  improvement  appear  to  follow  the  first  case. 

For  the  clear  day,  the  overall  effect  of  the  smaller 
regression  coefficients  was  to  reduce  the  estimated 
irradiance  values  in  the  model.  Figures  13-16  show  that 
the  revised  regression  coefficients  do  not  change  the  shape 
of  the  curve  for  the  calculated  values,  only  the  height  at 
which  the  curve  peaks.  This  helped  to  lower  the  morning 
and  evening  overestimations  but  it  also  added  to  the  mid 
day  underestimations.  This  resulted  in  higher  MBE%  and 
RMSE%  values.  The  magnitude  of  the  reduced  performance 
for  a  clear  day  is  small  and  since  the  majority  of  the 
days  in  the  study  are  cloudy  the  reduced  performance  of 
the  model  for  clear  days  is  acceptable. 

c.  Dav  232 

Nearly  every  site  experienced  one  or  two  days  that 
were  not  handled  well  by  the  model.  Satellite  analysis 
and  ground  observations  indicate  that  a  thin  band  of 
cirrus  traversed  across  the  state  during  day  232.  The 
daily  tables  indicate  that  day  232  was  a  problem  day 
shared  by  most  of  the  sites  in  the  network  and  that  the 
model  tends  to  considerably  overestimate  a  thin  cirrus 


condition . 
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The  daily  results  (Appendix)  for  day  232  indicate 
consistently  poor  MBE  and  RMSE  values  across  the  network. 
Correlation  coefficients  range  from  as  low  as  -.0515  at 
Logan  to  .8581  at  Springville  with  the  majority  of  the  R 
values  at  . 6  or  lower.  Due  to  the  consistently  poor 
results  over  the  entire  network,  it  appears  the  model  has 
some  problem  calculating  accurate  irradiance  values  for 
days  with  a  thin  cirrus  overcast. 

d.  Summary 

This  initial  test  of  the  Hay  and  Hanson  (1978)  model 
has  shown  that  the  regression  coefficients  developed  for 
the  tropical  Atlantic  were  not  completely  suitable  for 
the  data  used  in  this  study,  but  the  fit  is  remarkably 
good  considering  where  the  model  was  developed.  The  error 
is  due  to  the  bias  of  this  data  set  to  partly  cloudy  and 
overcast  conditions  and  the  higher  elevations  experienced 
over  the  network. 

Revised  regression  coefficients  lead  to  increases  in 
performance  for  partly  cloudy  sky  cover  conditions.  Al¬ 
though  overcast  conditions  were  never  experienced  over 
the  entire  network.  The  results  indicate  the  model  seems 
to  considerably  overestimate  a  heavy  cloud  condition. 

Under  clear  skies,  the  pattern  of  over  and  under  esti¬ 
mates  observed  in  the  initial  test  of  the  model  was  main¬ 
tained  with  increased  underestimation  during  the  mid  day 


hours.  This  leads  to  a  slight  decrease  in  performance  for 
clear  sky  conditions.  For  operational  use  of  this  model 
at  these  elevations  the  basic  shape  function  in  the  look¬ 
up  table  needs  to  be  modified. 

The  performance  of  the  model  for  partly  cloudy  condi¬ 
tions  showed  improvement  at  sites  that  were  cloudy  during 
the  mid  day  hours.  The  model  continues  to  experience 
trouble  calculating  irradiance  values  when  the  site  changes 
widely  from  hour  to  hour  due  to  cloud  shadowing  or  cloud 
movement . 

2 .  The  Tarpley  model 

An  aspect  of  some  importance  to  this  study  is  that  the 
model  be  applied  to  the  data  set  as  developed  by  the 
original  developer.  This  guarantees  that  any  differences 
in  results  cannot  be  attributed  to  changes  in  the  model 
and  must  reflect  the  different  environments.  There  are, 
however,  several  differences  from  the  original  Tarpley 
model  and  the  one  utilized  here. 

One  difference  between  the  two  models  is  the  brightness 
scale  used  in  the  satellite  data.  A  6-bit  scale  is  used  in 
the  Tarpley  (1979)  study  with  count  values  from  0-63.  In 
the  present  study  an  8-bit  scale  is  utilized  with  count 
values  from  0  -  255.  This  requires  that  several  values 
within  the  model  be  multiplied  by  a  factor  of  4  to  make  the 
two  scales  compatible. 
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An  additional  difference  is  the  source  of  the  precipi- 
table  water  values  necessary  for  the  irradiance  calcula¬ 
tions.  The  values  used  in  the  Tarpley  (1979)  study  were 
accessed  from  the  National  Meteorological  Center  (NMC) 

0000  GMT  operational  file.  Accessing  routines  retrieved 
precipitable  water  for  each  site  from  the  nearest  NMC  grid 
point  which  could  be  as  far  as  2  latitude  and  longitude 
from  the  station  (Tarpley,  1979)  .  In  this  study  the  pre¬ 
cipitable  water  values  were  calculated  for  each  hour  with 
Smith's  (1966)  empirical  formulation  using  temperature  and 
dew  point  data  available  from  each  site.  Although  the 
model  is  not  very  sensitive  to  precipitable  water,  this 
method  could  permit  a  more  accurate  assessment  of  the 
model . 

The  initial  test  of  the  Tarpley  (1979)  model  was  per¬ 
formed  over  all  the  test  sites.  Prior  to  the  test,  Julian 
days  172,  174,  and  233  were  selected  for  analysis  repre¬ 
senting  clear,  partly  cloudy  and  overcast  sky  conditions, 
respectively.  The  hourly  calculations  for  the  three  days 
were  compared  to  the  measured  irradiance  values  at  all 
the  sites  in  the  study.  Complete  statistics  for  the  sites 
are  listed  in  the  Appendix,  while  averaged  hourly  statistics 
for  a  few  select  locations  are  given  in  Table  7.  Figures 
24  -  33  depict  the  behavior  of  the  model  for  the  days  and 
sites  in  Table  7.  As  with  the  Hay  and  Hanson  (1978)  model, 
the  initial  test  indicated  no  bias  towards  any  of  the  test 
sites . 
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Table  7.  The  hourly  statistics  at  select  sites  using  the 
original  regression  coefficients  from  the  Tarpley  (1979) 


model 

233) 

for  clear  (day  172) 

sky  conditions. 

,  and 

partly 

cloudy 

(days  174, 

Day 

IRRADIANCE 

OBS  CALC 

N 

MBE% 

RMSE% 

LOCATION 

172 

26049.0 

26782.4 

11 

-2.7 

9.1 

CEDAR  CITY 

172 

25942.9 

26954.6 

11 

-3.8 

12.4 

ST.  GEORGE 

172 

26494 . 1 

26760.0 

11 

-1.0 

3.3 

GARLAND 

172 

26305.7 

26825 . 3 

11 

-1.9 

6.4 

SPRINGVILLE 

174 

20562.5 

22582.9 

11 

-8.9 

29.7 

CEDAR  CITY 

174 

23403.2 

23598.8 

11 

-.8 

2.7 

MILFORD 

174 

23564.4 

24777.9 

11 

-4.9 

16.2 

DELTA 

233 

12459.4 

14712.7 

9 

-15.3 

45.9 

MILFORD 

233 

15059.9 

20117.2 

9 

-25.1 

75.4 

PARK  CITY 

233 

14469.1 

18254 . 3 

9 

-20.7 

62.2 

LOGAN 

233 

19094.4 

21182.8 

9 

-9.9 

29.6 

GARLAND 

TOT: 

233404.6 

252548.9 

113 

-7.5 

The  averaged  hourly  calculations  from  the  initial  test 
indicate  that  the  model  generally  overestimates  the  observed 
irradiance  values  for  most  conditions. 


24.  Calculated  and  observed  irradiance  using  the  original  regre 
coefficients  from  the  Tarpley  (1979)  model. 


Calculated  and  observed  irradiance  using  the  original  regression 
coefficients  and  the  Tarpley  (1979)  model. 


LOCAL  TIME  (DST) 

r;  =  OBS  +  =  CALC 

Calculated  and  observed  irradiance  using  the  original  regression 
coefficients  from  the  Tarpley  (1979)  model. 


from  the  Tarpley  (1979)  model. 


Fig.  29.  Calculated  and  observed  irradiance  using  the  original  regre 
coefficients  from  the  Tarpley  (1979)  model. 


from  the  Tarpley  (1979)  model. 


from  the  Tarpley  (1979)  model. 


from  the  Tarpley  (1979)  model. 


from  the  Tarpley  (1979)  model. 


Calculated  and  observed  irradiance  using  the  original  regression 
coefficients  from  the  Tarpley  (1979)  model. 
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Figures  24-27  depict  the  trend  of  the  model  for  a  clear 
day.  About  sixty  percent  of  the  sites  show  curves  similar 
to  Figure  24.  The  figure  indicates  the  model  tends  to 
underestimate  the  irradiance  during  the  morning  and  mid 
day  hours  while  the  model  tends  to  overestimate  the  after¬ 
noon  and  evening  hours.  This  trend  would  produce  an  over¬ 
estimation  of  the  daily  irradiance  value. 

For  a  partly  cloudy  day.  Figures  28-34  indicate  that 
the  model  generally  follows  the  observed  trend  line  but 
usually  overestimates  the  early  morning  and  mid  day  condi¬ 
tions  and  often  misses  major  events  of  sunshine  or  cloudi¬ 
ness.  Again,  it  must  be  remembered  that  the  irradiance 
calculation  is  based  on  one  satellite  image  per  hour. 

This  may  not  represent  the  actual  condition  observed  at 
the  site. 

The  accuracy  of  the  original  model  was  determined  by 
running  a  statistical  analysis  of  data  listed  in  Table  7. 
The  totaled  value  listed  at  the  bottom  of  the  table  indi¬ 
cates  the  model's  performance  under  all  sky  cover  condi¬ 
tions.  Limitations  in  the  Tarpley  (1979)  model  become 
apparent  when  it  is  used  to  calculate  cloudy  conditions. 
Although  the  performance  of  the  model  is  quite  good  for 
clear  days,  an  attempt  is  made  to  improve  the  model's 
overall  performance  by  determining  new  regression  coef¬ 
ficients  . 


a.  Regression  coefficients 

The  development  of  new  regression  coefficients  was 
accomplished  using  a  subset  of  the  data  that  represented 
various  conditions  (elevation,  latitude,  longitude,  etc.) 
over  the  entire  network.  The  analysis  centered  on  days 
172,  174  and  233  to  represent  various  sky  cover  conditions. 

The  Tarpley  (1979)  model  contains  two  sets  of 
regression  coefficients.  The  first  set  of  values  correspond 
to  the  three  cloud  categories.  Since  the  original  model 
was  developed  for  a  similar  season,  elevation  and  predomi¬ 
nate  cloud  type  the  original  coefficients  are  utilized 
here.  The  other  set  of  regression  coefficients  are  related 
to  the  clear  brightness  calculation.  The  clear  brightness 
calculation  is  accomplished  for  each  hour  and  used  in  the 
final  irradiance  calculation.  Clear  brightness  (B)  is 
listed  as: 

2 

B  =  a  +  bcos  9  +  csm9  cost})  +  dsm9  cos  <t>  (Tarpley,  1979) 
where  9  is  the  local  solar  zenith  angle, <)>the  sun-satellite 
azimuth  angle  and  (a) ,  (b) ,  (c) ,  and  (d)  are  the  regression 

coefficients  of  interest. 

The  new  regression  coefficients  were  calculated  using 
an  available  software  package  called  Number  Cruncher. 
Coefficients  would  ideally  be  determined  for  each  study 
site,  although  Raphael  and  Hay  (1984)  showed  the  use  of 
locally-derived  coefficients  did  little  to  improve  the 


performance  of  the  model.  The  revised  coefficients  are: 
a  -  45.32,  b  =  61.74,  c  =  12.85,  d  =  22.12. 

b.  The  revised  model  results 

The  revised  model  was  tested  over  all  the  sites  using 
all  the  available  data.  Statistical  comparison  is  conducted 
on  the  same  days  used  in  the  initial  test  of  the  model. 

The  hourly  average  statistics  at  select  sites  are  listed 
in  Table  8  (see  Appendix  for  statistics  on  the  individual 
sites)  . 

From  the  results  in  Table  8  and  the  Appendix,  an 
improvement  in  the  model  was  observed  for  partly  cloudy 
sky  conditions.  As  with  the  Hay  and  Hanson  (1978)  model, 
a  slight  decrease  in  performance  is  seen  for  cloudy  days, 
a  result  of  the  bias  of  our  data  toward  cloudy  days. 

Figures  35-38  indicate  the  performance  of  the  model 
for  a  clear  day  using  the  revised  regression  coefficients. 
The  morning  and  mid  day  estimates  are  now  very  close  to  the 
observed  irradiance  values  but  the  modelte  performance  tends 
to  slightly  decrease  during  the  afternoon  and  evening  hours. 
This  suggests  that  this  trend  may  be  due  to  an  increasing 
solar  zenith  angle  and  a  high  aerosol  count  that  would  be 
be  associated  with  a  clear  afternoon  rather  than  a  clear 
morning . 

Figures  39-45  indicate  that  the  revision  of  the 
regression  coefficients  can  effect  the-  shape  of  the 


Table  8 


The  hourly  statistics  at  select  sites  using  the 


revised  coefficients  from  the  Tarpley  (1979)  model  for 
clear  (day  172),  and  partly  cloudy  (days  174  and  233)  sky 
conditions . 


DAY 

IRRADIANCE 

OBS  CALC 

N 

MBE% 

RMSE% 

LOCATION 

172 

26049.0 

25441.3 

11 

2.4 

7.9 

CEDAR  CITY 

172 

25942.9 

26552 . 5 

11 

-2.3 

7.6 

ST.  GEORGE 

172 

26494 . 1 

26325.1 

11 

.6 

2.1 

GARLAND 

172 

26305.7 

26394.6 

11 

-.3 

1.1 

SPRINGVILLE 

174 

20562.5 

19351.4 

11 

6.3 

20.8 

CEDAR  CITY 

174 

23403.2 

18830.6 

11 

24.3 

80.5 

MILFORD 

174 

23564.4 

22393.8 

11 

5.2 

17.3 

DELTA 

233 

12459.4 

13078.2 

9 

-4.7 

14.2 

MILFORD 

233 

15059.9 

18697.9 

9 

-19.5 

58.4 

PARK  CITY 

233 

14469.1 

16802.6 

9 

-13.9 

41.7 

LOGAN 

233 

19094.4 

19954.4 

9 

-4.3 

12.9 

GARLAND 

TOT: 

233404.6 

233822.4 

113 

-.2 

estimated  irradiance  value  curve,  rot  just  the  height  at 
which  the  curve  peaks.  The  calculated  line  now  mirrors  the 
observed  irradiance  line  quite  well.  The  model,  however, 
has  a  tendency  to  underestimate  the  observed  irradiance 
value.  Occasionally,  the  performance  of  the  model  also 


Calculated  and  observed  irradiance  using  the  improved  regre 
coefficients  in  the  Tarpley  (1979)  model. 


DAY  172  C 


using  the  improved  regression 


Fig.  38.  Calculated  and  observed  irradiance  using  the  improved  regression 
coefficients  in  the  Tarpley  (1979)  model. 
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Fig.  41.  Calculated  and  observed  irradiance  using  the  improved  regression 
coefficients  in  the  Tarpley  (1979)  model. 


i 

9 


(unoH/t-fi  nO  jawvwvw 


Fig.  44.  Calculated  and  observed  irradiance  using  the  improved  regression 
coefficients  in  the  Tarpley  (1979)  model. 
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Fig.  45.  Calculated  and  observed  irradiance  using  the  improved  regression 
coefficients  in  the  Tarpley  (1979)  model. 
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appears  to  suffer  from  the  inability  of  only  one  satellite 
image  to  represent  the  conditions  over  the  entire  hour. 
Additionally,  the  model  continues  to  have  some  trouble 
with  a  thin  overcast  day  (day  232),  as  with  the  Hay  and 
Hanson  (1978)  model. 

3 .  Summary 

The  summary  statistics  for  clear,  partly  cloudy,  and 
overcast  sky  conditions  and  the  total  sample  are  presented 
in  Table  9  for  both  of  the  revised  models.  Several  factors 
should  be  noted. 

The  results  from  this  study  indicate  that  the  Hay  and 
Hanson  (1978)  model  performs  slightly  better  than  the 
Tarpley  (1979)  model  for  both  the  total  hourly  and  daily 
values.  Raphael  and  Hay  (1984),  reported  similar  results 
from  their  use  of  the  two  models  over  southwest  Canada. 

The  Tarpley  (1979)  model  performed  better  for  a  clear 
day  while  the  Hay  and  Hanson  (1978)  model  performed  better 
for  the  partly  cloudy  situation. 

For  the  Hay  and  Hanson  (1978)  model  the  daily  (RMSE% 
9.7)  and  hourly  (RMSE%  19.1)  statistics  for  the  total 
sample  were  similar  to  those  quoted  by  Hay  and  Hanson 
(1978)  . 

For  the  Tarpley  (1979)  model,  the  daily  correlation 
coefficient  of  .9611  was  slightly  better  than  the  one 
quoted  by  Tarpley  (1979).  The  short  term  performance 


Table  9.  Summary  statistics  showing  the  performance  of 
the  original  and  revised  models  for  clear,  partly  cloudy 
and  overcast  sky  conditions  on  an  hourly  (h)  and  daily  (d) 
basis.  Units  are  KJ  m-2/hour. 


Condition 

N 

Time 

RMSE 

RMSE 

(%) 

Mean 

Obs  Calc 

Revised  Hay  Model 

Clear  days 

294 

h 

234.2 

9.4 

2470.2 

2438.9 

Partly  Cloudy 

519 

h 

421.5 

18.3 

2221.8 

2302.9 

days 

Total  sample 

1183 

h 

411.7 

19.1 

2043.4 

2159.3 

137 

d 

1783.5 

9.7 

17371.0 

18302.3 

Revised  Tarpley  Model : 

Clear  days  293 

h 

240.8 

9.4 

2482.1 

2542.4 

Partly  Cloudy 

508 

h 

516.6 

21.8 

2223.3 

2367.4 

days 

Total  sample 

1151 

h 

518.9 

23.9 

2060.3 

2172.7 

137 

d 

1949.8 

10.7 

17219.9 

18150.8 

statistic  obtained  for  the  total  sample  is,  however,  higher 
than  that  reported  by  Tarpley  (1979)  for  both  the  hourly 
and  daily  analysis. 


CHAPTER  V 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  intention  of  this  thesis  was  to  determine  if  a 
reliable  method  of  estimating  solar  irradiance  could  be 
demonstrated  for  a  mountainous  region  using  satellite 
data.  A  comparison  has  been  accomplished  using  two  statis¬ 
tical  models,  the  results  showing  slightly  better  perfor¬ 
mance  from  one  of  these  models.  Some  anomalies  were 
expected  and  found.  The  development  of  regression  coef¬ 
ficients  more  suitable  to  this  region  improved  the  perfor¬ 
mance  of  both  models  for  most  situations.  In  comparing 
the  results  of  this  study  with  the  work  done  by  the  original 
developer,  similar  results  were  obtained. 

The  results  from  this  study  demonstrates  slightly 
better  performance  by  the  Hay  and  Hanson  (1978)  model 
over  the  Tarpley  (1979)  model  in  estimating  solar  irradiance 
under  partly  cloudy  and  overcast  skies.  The  Tarpley  (1979) 
model  proved  better  able  to  estimate  irradiance  under  a 
clear  sky.  Both  models  systematically  overestimate  over¬ 
cast  conditions  while  for  partly  cloudy  conditions,  the 
Hay  and  Hanson  (1978)  model  generally  underestimates  and 
the  Tarpley  (1979)  model  overestimates. 

Several  irregularities  were  obesrved  at  the  Park  City 
site,  the  only  site  that  is  truly  in  mountainous  terrain. 
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The  effects  of  high  elevation  and  complex  terrain,  no  doubt 
cause  variability  in  both  models.  Sites  at  higher  eleva¬ 
tions  characteristically  receive  more  insolation  than  sites 
at  low  elevation.  The  development  of  new  regression  coef¬ 
ficients  may  help  to  alleviate  some  of  this  variability. 

The  simple  statistical  framework  of  both  models  makes  their 
own  use  over  very  complex  terrain  a  difficult  matter. 

As  pointed  out  by  the  original  developers,  the  models 
have  several  shortcomings.  Both  models  inadequately  handle 
cloud  absorption.  It  has  been  estimated  that  up  to  25%  of 
the  incident  visible  radiation  can  be  absorbed  by  certain 
clouds  (Tarpley,  1979) .  It  would  be  necessary  to  determine 
cloud  type  and  thickness  to  account  for  this  factor. 

Another  factor  to  be  considered  for  inclusion  within  the 
models  are  the  effects  of  aerosols. 

For  particular  use  over  a  high  elevation  region  with 
geographic  variations  in  albedo  -  such  as  Utah,  several 
recommendations  are  in  order.  The  look-up  table  calibration 
technique  used  in  the  Hay  and  Hanson  (1978)  model  needs 
further  study.  For  the  Tarpley  (1979)  model,  the  regression 
coefficients  used  in  determining  the  three  cloud  categories 
may  not  be  totally  suitable  for  this  area.  The  precipitable 
water  values  calculated  for  the  study  did  not  take  into 
account  the  high  elevations  experienced  over  the  network. 

The  precipitable  water  values,  therefore,  may  be  too  large. 
The  results  produced  from  both  models  depend,  to  a  great 
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The  following 
Appendix . 

OBS  = 
CALC  = 

N  = 
MBE  = 
RMSE  = 
MBE%  = 
RMSE%  = 
R  = 
RSQ  = 
JD  = 


notation  applies  to  the  Tables  listed  in  the 


Observed  irradiance  summed  for  the  day. 

Calculated  irradiance  summed  for  the  day. 

Number  of  hours  used  in  the  calculation. 

Mean  bias  error 

Root  mean  square  error 

Relative  mean  bias  error 

Relative  root  mean  square  error 

Correlation  coefficient 

Coefficient  of  determination 

Julian  day 


19432.4  18427.1  8  103.2 


17292.4 
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